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1. Introduction 

Monte Carlo methods have tremendously improved our understanding of a variety of strongly 
interacting quantum field theories. Beginning with the discovery of "cluster algorithms" for clas- 
sical spin systems ^ and the "loop algorithms" for quantum spin systems ^ the search for 
efficient algorithms has been an important area of research over the past two decades. Although 
many algorithmic improvements have been achieved for QCD-like problems [Q], it has been espe- 
cially difficult to find good algorithms for lattice field theory problems in the presence of a chemical 
potential or when the microscopic theory contains strongly interacting massless fermions. How- 
ever, recently the landscape has begun to change somewhat. The progress over the past decade, 
especially in studies of quantum spin systems and strong coupling lattice gauge theories, suggests 
the existence of a new approach to a class of lattice field theories. In this talk I will outline the 
ideas behind the new approach by discussing a few examples. 

One essential feature of the new approach is to formulate lattice field theories in representa- 
tions that are similar to the world-line representation. Such representations arise naturally in the 
Hamiltonian formulation but are not usually used in the Lagrangian formulation. As we will see 
below, they also arise naturally in Euclidean lattice field theory when one thinks about the strong 
coupling (or high temperature) expansions. In certain cases the representations are indeed the well- 
known world-line representations, in others they look different but share many properties of the 
world line representations. For this reason we will refer to the new approach generically as the 
"WL-approach". As we will see below one of the advantages of the world-line representations is 
that there are no new sign problems because of adding a chemical potential unlike the conventional 
approach. On the other hand, the world-line representations are not friendly for local Monte Carlo 
methods due to the presence of constraints. The loop cluster algorithm was the first successful 
algorithm that was able solve the constrained problem efficiently. However, this algorithm was 
rather restrictive and it was difficult to maintain the efficiency in the presence of different types of 
couplings especially with a chemical potential. Recently, it was realized that the loop algorithm 
is just one example of a more general class of algorithms refened to as "worm algorithms" [||] or 
"directed-loop algorithms" [§] or "directed-path algorithm" [0]. We will refer to these slightly dif- 
ferent algorithms together as "worm-algorithms" based on the pioneering work of the first authors 
Prokofev and Svustinov. These algorithms introduce a defect in the constraints and sample the de- 
fect space until the defect naturally eliminates itself. The sampled defect space is itself a measure 
of the two point correlation function. When combined with the worm-algorithm, the WL-approach 
to lattice field theories have been found to be as efficient as the original cluster algorithms for clas- 
sical spin models [Q, |8||. On the other hand the WL-approach seem to be more widely applicable. 
As we will discuss below, the worm-algorithm can in principle be applied to CP^ ^ models which 
cannot be solved using the conventional cluster algorithms. 

As with any Monte-Carlo method, the WL-approach can fail due to sign problems. Sign 
problems can arise due to frustrating interactions or in the presence of fermions. In two dimensions 
fermionic sign problems can sometimes be solved in the world line representation since fermions 
resemble hard-core bosons [^, |l^]. Thus, WL-approach should be applicable to these problems. 
Solutions to fermionic sign problems also emerge in strongly coupled lattice gauge theories with 
staggered fermions. In this case configurations naturally resemble world-line representations where 
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fermions are confined into bosons. Wlien combined with the efficient directed-path algorithm 
this provides a new and powerful WL-approach for some scalar field theories which include those 
that are invariant under symmetries like SU{l)xU{\) and SU {2) x SU{2) xU{l) ^^,^]. 
These theories, which are of interest in the context of understanding QCD-like theories, could not 
be solved efficiently with conventional approaches. 

In higher dimensions solving fermionic sign problems is difficult, however the WL-approach 
suggests new solutions in many cases [11, 12, [13, 14 1. If the sign problems can be solved in the 
WL-approach, then one can combine worm-algorithms with other conventional algorithms to solve 
the problem. Below we will discuss one such algorithm for the massless lattice Thirring model in 
any dimension. In the WL-approach the partition function is rewritten such that it resembles the 
well known bag-model of QCD where the bags are dynamically determined objects. Inside the bag 
the fermions are free, while outside they are confined. For this reason we call this new fermion 
algorithm the "dynamical-bag algorithm". 

The basic ideas of the WL-approach are also applicable to gauge theories. Here one needs to 
rewrite the problem in a representation of world-sheets giving us a new "WS-approach" where a 
worm-type algorithm involves cutting and reconnecting world sheets. Although very little is known 
about these surface-cluster algorithms, the ideas seem promising. Below we will show some results 
which suggest that it may be possible to measure large Wilson loops with little effort in the confined 
phase in Abelian lattice gauge theory using a WS-approach. 



2. World-Line Approach: XY Model with a Chemical Potential 

In this section we will illustrate the WL-approach using the example of the XY model. To 
make things interesting we will introduce a chemical potential that couples to the U{1) charge of 
the particles. We begin with the action of the lattice field theory, in the conventional formulation, 
which is given by 




(2.1) 



where x is a point on a J + 1 -dimensional hypercubic lattice, a = t,l,2..,d represents the direction, 
j8 is the inverse couphng of the theory and exp{i^x) is the non-linear U{1) bosonic field. Clearly, 
when jJ. ^0 the action is complex and the path integral suffers from a sign problem and one cannot 
use the conventional Monte Carlo methods. In particular the Wolff-cluster algorithm ^ is no 
longer applicable. 

In order to construct the WL-approach we rewrite the partition function using the strong cou- 
pling (or high temperature) expansion to all orders. Using the identity 

exp{j8cos(0)}= £ 4(J8)e'*^^ (2.2) 

on every bond and integrating over the original angle variables the partition function can be written 
in terms of integer bond variables kx,a and one gets 
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Figure 1: An illustration of the world-line configuration for the XY-model. 



Here 4(j3) is the modified Bessel function of the first kind. The delta function in the above expres- 
sion shows that the bond variables satisfy a local constraint at each lattice site x which is nothing 
but the current conservation relation. Thus, the partition function of the non-Unear sigma model 
has been rewritten as a sum over configurations of current-loops. Figure ^ illustrates one such con- 
figuration in two dimensions. Importantly, this representation does not suffer from a sign problem 
even at non-zero /i. The sign problem has been traded for a constraint condition which is not easily 
satisfied by a local change. However, developments over the past decade have shown that there are 
efficient non-local updates for such constrained problems ^, ^ [T8| ]. 

Is this approach general? In other words can the partition function of a lattice field theory be 
written in a world-line representation and then a worm algorithm be used to solve the problem? We 
will argue below that the answer is "yes" in many cases. This was recently shown to be true for the 
Ising model [^. In the next section we will argue that the CP^^^ model can also be solved in this 
approach. 



3. World-line Approach: CP^^ Model 

In this section we will construct the world-line representation of the CP^ ^ model as a non- 
trivial application of the WL-approach. One may recall that conventional cluster algorithms do not 
work for these models even in the absence of a chemical potential. CP^ ^ models were recently 
formulated as the low energy effective theory of an SU (A/^) -symmetric quantum-spin system in 
the Hamiltonian (D-theory) formulation [p^. The spin system was formulated in the world-line 
formulation, and solved using the loop-cluster algorithm. Thus, the first successful efficient cluster 
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algorithm for CP'^^ ^ models emerged in the WL-approach. Here we will argue that a world-line 
representation can also be constructed for the conventional model directly in Euclidean space. 

The CP^^^ model is conventionally written in terms of A'^ component complex vector zf,a = 
1,2, ..,N at each lattice site / with the constraint |z?p = 1 The action of the model is 

given by 

5=-^ (3.1) 

<'j> 

where < ij > denotes the bond connecting the nearest neighbor sites / and j. The action is invariant 
under global SU{N) transformations Zi Uzi and U{1) gauge transformations Zi exp(/0,)z,-. 
Since the gauge field is not dynamical in the theory, the model is equivalent to strongly coupled 
A^-flavor scalar QED. Let us now derive the world-line representation. 
Consider the partition function of the model, 

Z= [ Ui^ZijcMP I {^zr)(j'jz"j)) (3.2) 

■■' i <ij> 

where [dzi] is the integral over 2A'^-dimensional unit vectors. The partition function can be rewritten 

as 

z = Jm-^] n n I (3-3) 

where each bond < ij > contains anN xN matrix of non-negative integers riij = n"j which deter- 
mines the power of (z^zf ) and {z'jZ")- Now one can integrate over [dzi] using the identity 

27i^ki\k2\h\...kN^. 



+ k2 + ... + kN+N-l)\ 

If we define k = {ki,k2,...,kf^) as an A/^- vector and 



(3.4) 



^ 27l''ki\k2\k,\...kN\ 

'^"^ {k,+k2 + ... + ki, + N-l)\' 
the partition function can be written compactly as 

^=ifnw)5,,inn£f) (3-6) 

[nij] V ' J \<U> ab ^ij- J 

where 

^i=IlK+.)+«(t,).}' /'i = Il{"S^+.)+<.)] (3-7) 
M * Pi b ^ ' 

Thus, the partition function has been written in terms of the constrained bond matrices nij. These 
constrains encode the conserved currents of the CP^^ model and hence we call this the world-line 
representation. Using the ideas developed in it should now be possible to update the constrained 
system. This has not yet been accomplished but is an interesting research project for the future. 
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Figure 2: An illustration of a closed pack dimer configuration in two dimension. In d+ 1 -dimensions 
statistical mechanics of such configurations naturally lead to the physics of the (i-dimensional XY model. 



4. Bosons as Fermionic Composites 

Scalar field theories are naturally formulated using scalar fields. However, all scalar particles 
discovered in nature until now are bound fermionic composites. Interestingly, we have discovered 
recently that it is also computationally simpler to formulate certain scalar lattice field theories from 
a microscopic model of fermions such that the bosons arise as fermionic composites. Composite 
fermion models of scalar field theories can be constructed easily as strongly coupled lattice gauge 
theories. They lead to a novel world-line approach where the fermionic nature restricts the number 
of allowed scalar particles on a given lattice site. The underlying idea is very similar to the D- 



theory formulations of field theories []22|], however the formulations discussed below occur directly 
in Euclidean space and often it is easier to design worm-type algorithms in this approach. Let us 
illustrate this with two examples. 

First consider the XY model. A simple XY model of composite fermions on a J + 1 dimen- 
sional hypercubic lattice is given by 



(4.1) 



xj=l,2-,d 



where ¥x,¥x ^'^'^ Grassmann valued lattice fields. This theory has an exact global U{1) sym- 
metry where Yx —>■ e'^^'^^Yx and \j7^ —>■ e"^"^!}/^ where Ox is +1 on even sites and —1 on odd sites. 
Note that one of the dimensions (referred to as t here) has been singled out and acts like a fictitious 
temperature. If we study this theory on a L"^ x Lf lattice with Lf fixed, then by tuning T one obtains 
the physics of the conventional XY model in d dimensions through dimensional reduction. 
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It is possible to integrate out the Grassmann variables and write the partition function as a sum 
over all possible closed packed dimer configurations. In particular every dimer in the temporal 
direction carries a weight T . Thus, the statistical mechanics of closed packed dimers on a cubic 
lattice in <i + 1 -dimensions leads naturally to the physics of the XY model. A configuration of a 
closed pack dimers is shown in Fig. |^ Clearly these configurations are much simpler to represent 
on a computer than the world-line configurations of the conventional XY model shown in Fig. |l|. 
Another advantage of this representation of the XY model is that it can be updated very efficiently 
with a worm-type algorithm which is much simpler to code than the Wolff-cluster algorithm 
It has been shown that this model and its variants do indeed reproduce the expected KT physics 



in J = 2 [23] and XY universal critical behavior in J = 3 [24]. Interestingly, in 1 + 1 dimensions 
with r = 1 and Lf = L, the model leads to a non-compact version of the XY model and so contains 
no vortices and evades a KT transition [^]. Thus, the fermionic representation also allows us to 
formulate naturally a non-compact U{\) field theory. 

Second we focus on an SU{2) x SU{2) xU{l) model of composite fermions. This is an 
interesting toy model for pions of two flavor QCD. The action of the model is given by 

S=- ^ Tr [I, !,+,•] -r^Tr[£,r,+,]-c£detZ;, (4.2) 

x,i=l,2..,d X X 

where 

Sv^l 7 ](ux d.] = \ Z- T-: I, (4-3) 




is a 2 X 2 matrix made up fermionic bilinears. Here Ux,Ux,dx and dx are four independent Grass- 
mann variables defined on each site. When c = 0, the action is symmetric under the SU (2) x 
SUS{2) xU{\) chiral transformations 

LLxR^d^ for X even (4.4) 

ELxL^Q-'^ for X odd 



When c / 0, the symmetry under the U{\) subgroup is broken explicitly. Hence the parameter c 
behaves like the anomaly. 

Once the Grassmann integration is performed the partition function can be written as a sta- 
tistical mechanics of two different types of closed oriented loops which form the world lines of 
the four pions in the model. A non-zero c creates lattice sites which the pion loops cannot touch. 
An example of a configuration of pions is given in Fig.|3[ We refer to such sites as "instantons". 
The model has been studied recently with a directed path algorithm and we refer the reader to the 
original work [[T^]. Cluster algorithms in the conventional formulation do not work for scalar field 
theory with an SU (2) x SU {2)xU{V) field theory. For this reason until recently it was difficult to 
study the effects of the anomaly on the two-flavor chiral transition. Using the worm-approach we 
were able to study this question non-perturbatively and beyond mean field theory for the first time 
and found that a rather strong anomaly was necessary before the 0(4) universality sets in [p^]. 
Using this model we have also clearly shown that the cr-resonance can have a significant impact on 
the width of the region where chiral perturbation theory is valid [26]. A variant of the above model 
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Figure 3: An illustration of a world-line configuration for an SU (2) x SU (2) x f/ ( 1 ) scalar field theory that 
arises from a composite fermion model. 

also allows one to study an SU (2) x ?7(1) scalar field theory which naturally arises from two-color 
QCD with staggered fermions [|T5]]. 

Lattice QCD with staggered fermions at a finite baryon chemical potential in the strong cou- 
pling limit can also be formulated and studied in the WL-approach. This was already done almost 
two decades ago [|^]. The theory however suffers from a sign problem which can only be solved 
at zero chemical potential. In any case, in the previous study the problem was solved using a brute 
force technique but using a local algorithm. Recently, the problem has been revisited but now us- 
ing the worm algorithm. Clear advantages of the worm approach has been demonstrated. For more 
details we refer to the work presented by Michael Fromm at this conference. 

5. World-line Approach: The Thirring Model 

Let us now discuss how the WL-approach can be applied to a fermionic field theory if the 
fermion sign problem can be solved. Solution to the fermion sign problem in three or more di- 
mensions usually involves a re-summation over a class of configurations. This makes the problem 
computationally more demanding than bosonic field theories since a significant fraction of the 
update-time goes into the re-summation effort. Clearly, the new approach does not alleviate this 
problem. On the other hand it offers new techniques to solve the sign problem and some of them 
are better than the conventional methods. Conventionally, the solution to the sign problem in the 
presence of fermions comes from computing the determinant of a matrix whose size grows with 
system size. On the other hand in the WL-approach other interesting solutions which are compu- 
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tationally less demanding can emerge. One such approach called the meron-cluster approach was 
found about a decade ago [^]. In the following we will device a new determinantal solution to the 
lattice thirring model. The main advantage of this solution is that the size of the matrix is dynami- 
cally determined by the parameters of the model and does not always grow with system size. This 
is particularly true deep in the strong coupling phase (confined phase) which is the phase where the 
conventional algorithms become more demanding. 

The action of the model we consider is constructed with massless staggered fermions and is 
given by 

S = -Y.^ii {^)'WxWx+ii - Wx-ii] - UWxVxWx+fi Vx+fi ■ (5. 1) 

Here we assume x is a point on an L'^ hypercubic lattice with toroidal boundary conditions, the index 
/I = 1,2,.., J labels the direction, the fermion fields Yx and \j7^ are single component Grassmann 
fields with anti-periodic boundary conditions. Since the interactions exist on bonds the model is 
equivalent to the massless Thirring model in the weak coupling limit [|^, 28]. The lattice model is 
invariant under an exact U{1) chiral symmetry: 

Yx exp(/a^0)i/A^.; \j7^ exp{-iax(l))Yx (5-2) 

where CJ^ is 1 on even sites and —1 on odd sites. For small U this chiral symmetry remains unbro- 
ken and the system contains massless fermions, while for large U the symmetry is spontaneously 
broken and the system contains massless Goldstone bosons but massive fermions. Thus, there is an 
interesting quantum phase transition in the thermodynamic limit at a critical coupling Uc- A variant 



of this phase transition may be of interest in graphene-like systems [ ]29| ] 



The conventional Monte-Carlo approach is to rewrite the partition function 

[d\j7dY]e)<jp{-S) (5.3) 



by introducing new auxiliary variables on the links connecting x with x + jj. so that the four- 

fermion operator is converted to a fermion bilinear. Mathematically one can show that if we begin 
with 

Z = j d4 [dYdY] exp| £"iFx(^[0]).v,yV^y} (5-4) 



x,y 



where M((/)) is an x matrix given by 



(5.5) 



and integrate over the angles < (p^ (x) < 2% one can recover the original partition function given 



in eq.(5.3). Bit if we integrate over the Grassmann variables first we obtain 

Z = j [d<p] Det(M(0)) (5.6) 

where the determinant is non-negative and thus the sign problem is solved. One can now use either 
the determinantal Monte-Carlo algorithm or the Hybrid Molecular Dynamics algorithm to generate 



[0] accordingly |27, 28 1. In the latter approach it is common to introduce a fermion mass in order 
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Figure 4: An illustration of a "dynamical-bag" configuration as discussed in the text. 

to regulate the condition number of the fermion matrix M[0]. However, this breaks the chiral 
symmetry and it is difficuh to extrapolate to the massless limit. The exactly massless problem of 
interest is known to be computationally very demanding especially for large values of U due to the 
excessively large number of small eigenvalues of M[<^] [|2^]. 

The above approach is the most commonly used method to deal with four-fermion couplings 
and is known in the literature as a Hubbard-Stratanovich transformation. Let us now discuss a WL- 
approach to the problem. Instead of introducing auxiliary variables, we begin with the partition 
function given by 

Z = j [d\\fd-W]Q^V [^ti{x)'WxWx+n - Wx-i^+UWxWxWx+^Wx+n]^ , (5.7) 
and expand it in powers of U using 

^^ViUWxWxWx+llWx+ll) = 1 + UW:,Vx'Wx+n¥x+ll- (5.8) 

The Grassmann integration then gives 

^= E (llU"''W<W[n]) (5.9) 

where = 0, 1. The Ux^^ = 1 bonds are the same as the hard core dimers encountered in the 
previous section in the fermionic XY model. Note that in this approach the Grassmann integration 
leads to a determinant of a different matrix W[n\, which is just the free fermion matrix where the 
sites connected to Ux^^ = 1 are dropped. It is easy to argue that Det(lV[?i]) is also non-negative and 
thus there is again no sign problem. 
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Figure 5: Plot of the number of n,,^ = 1 bonds on a 16^ lattice at C/ = 1 in the dynamical-bag algorithm. 



Interestingly, the configuration [n] divides the lattice into "bags" of sites connected with only 
nx.n = bonds. Inside a bag the fermions hop freely while outside they are confined to hop on 
single bonds. The size and shape of the bags are dynamically determined by the value of U . Figure 
^ gives an illustration of one such configuration. For this reason we call our new algorithm the 
"dynamical-bag" algorithm. Since a single world line configuration of fermions inside the bag can 
give negative weights due to the Pauli principle, one has to resum all the possible fermion world 
lines within the bag. This gives the Det(W[?i]), which in our case is non-negative. When U is large, 
the bags are small comprising of a few neighboring bonds and thus independent of the system size. 
Hence, the computation of the determinant must be easy. In fact, as we will argue below, it is 
not even necessary to compute the determinant, but just one matrix element of the inverse of 
at every local update step, which is quite easy for small matrices. Note that it is precisely when 
U is large that conventional algorithms begin to fail. When U is small the bags can percolate and 
become as big as the system size. This naturally leads to massless fermions. Here the WL-approach 
will be similar to the conventional approach in efficiency. The typical size of the bags may provide 
an interesting length scale for the problem. 

Our algorithm consists of two steps: one step changes Wj^^ between and 1 on the bonds and 
the other step moves the bonds n^,^ bonds around. One can develop a combination of a regular 
local heat-bath algorithm to accomplish the former and a worm-type algorithm to accomplish the 
latter. An important point to note is that the probability to introduce or remove a bond depends 
on the ratio of two determinants one with and one without the bond. This ratio is exactly equal to 
|[(W[?i])^'].:c,.v+/jP where the bond in consideration is between the sites x and x + /i and the matrix 
W[n] does not contain the bond under consideration. The computation of the inverse can be a 
time consuming process. Further, if the matrix W[n] contains exact zero modes, one has to know 
about its existence so that such configurations are not generated. These two difficulties makes the 
algorithm more time consuming as compared to bosonic algorithms. However, this is a price one 
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has to pay for doing fermionic physics and we do not know a way out of it at the moment. However, 
at large U when the "fermions" are mostly bound into bosons these difficulties become milder and 
disappear completely at infinite U . Thus, the algorithm knows when to work hard to take into 
account the dynamics of the fermions. In Fig. || we show the time evolution of the number of 
nx,iL = 1 bonds on the lattice at ?7 = 1 starting from two different initial conditions on a 16^ lattice 
which was accomplished in a few hours on a laptop with just a local heat bath algorithm. 

Although the local heat bath algorithm may be sufficient at small U , it will become inefficient 
as U increases. At large U the worm part of the algorithm is essential for efficiency. The worm al- 
gorithm is also necessary to measure observables such as the chiral condensate susceptibility which 
can get contributions from fermionic correlations between two different bags. These correlations 
are difficult to measure in the WL-approach using simple matrix inversions. In the worm algorithm 
such correlations come from configurations which contain two defects (the head and the tail of 
the worm) one of which is present in one bag and the other in the other bag. We are currently 
implementing this complete algorithm and the technical details will be discussed elsewhere. 



6. World- Sheet Approach: Abelian Gauge Theory 

The natural extension of world-lines to gauge theories must be in the form of world-sheets. 
Thus, the WL-approach must be modified into a WS -approach when applying to gauge theories. 
Here we illustrate the ideas using a pure compact U{\) lattice gauge theory in four dimensions. 
In the conventional formulation the theory describes the dynamics of the compact gauge field < 
^a{x) < 2n,a = 1,2,3,4 which are angular variables associated to bonds {x,a) that connect sites 
X andx + a. The model is described by the partition function 



[#]exp[-S]. (6.1) 
where the action is taken for convenience to be the Wilson action 

S=-I5 £ cosf^^vW), (6.2) 

where the sum is over all plaquettes P{x,}JL, v) (defined by a lattice site x and two forward directions 
/I < v) and 0^vW = 0^1 W + <^v{x + ci(l) — <^p,{x + av) — (j)v{x). We define (j)-a{x) = —<^a{x — a) 
for convenience. 

Motivated by the world-line approach we can rewrite the partition function in terms of world- 
sheet variables. For this purpose we begin with 

Z= j [d(^] Y\ e/^™^WMvW). (6.3) 



and use eq.(2.2) on every plaquette and then integrate over all angles 0a(x) as in the XY model. 
We now get 

2 = L n 4^vW(^) n ^v[^MvW-'tMv(^-v)],0 (6-4) 

[k\ P{x,ii,v) {x,ti) 
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where [k] represents a configuration of constrained integers on each plaquette /"(x,/!, v). For sim- 
plicity we define k^y{x) = —ky^{x). The configurations must now satisfy the constraint that 

£[/:^v(x)-fc^,v(x-v)]=0 (6.5) 

V 

at every bond (x,/x) and are shown with delta functions in the above expression. It is possible to 
argue that such a constrained set of integers arise from oriented closed surfaces or world sheets. 
For example, an update of [k\ requires the change of ^ on a closed surface of plaquettes. Any 
configuration [k\ can be constructed by a series of local cubical updates and sheet updates in two 
dimensional planes. 

An interesting observable is the Wilson loop. Consider a closed loop on a /XV plane of size 
Lx L. The Wilson loop associated with this loop is defined as 

N 

WL = llexpii^a,ixk)) (6.6) 

k=i 

where Xk and at are the sites and directions along the loop. The average of this observable is given 
by 

where the final average is assumed to be performed in the world-sheet representation and P stands 
for the set of plaquettes that cover the surface of the Wilson loop. 

In order to explore the usefulness of the new representation, we have implemented a worm- 
like update for the [k] configuration which involves picking a /iv plane at random and introducing 
a defect in the world-sheet configuration in the form of a 1 x 1 Wilson-loop. We then increase or 
decrease the size of the defect-loop but restricting it to be a square loop within the plane. During this 
process we allow for local cubical updates within a slice on either side of the sheet, which allows 
the surface being updated to fluctuate. The algorithm ends when the defect closes on itself or 
propagates and updates a whole sheet. We have observed that sheets do get updated in the Coulomb 
phase and when the volumes are small. Figure ^ shows a comparison of the results of the world- 
sheet algorithm with the conventional algorithm for various sizes of a square shaped Wilson-loop. 
The figures show that for small sizes both algorithms are comparable and produce the expected 
physics, namely an area law in the confined phase and a perimeter law in the Coulomb phase. 
Interestingly, the world-sheet algorithm is able to measure large Wilson-loops quite well in the 
confined phase, while the conventional algorithm still appears superior in the Coulomb phase. 

7. Conclusions 

We have tried to demonstrate in this talk that there is a new computational approach to a variety 
of lattice field theories. This approach is based on rewriting the partition function in a world-line 
representation. The constrained configurations in this representation can be updated efficiently with 
the newly discovered worm-algorithms. This approach gives an alternative method for studying 
many scalar field theories. They also lead to new solutions to the fermion sign problems. One such 
solution leads to an interesting fermion algorithm which we call the dynamical-bag algorithm. 
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Figure 6: Comparing the results for the square Wilson-loop as a function of its size obtained using the 
world-sheet algorithm (filled circles) and the conventional algorithm (open squares). The world-sheet algo- 
rithm is better than the conventional method in the confined phase (left figure) while the opposite is true in 
the Coulomb phase (right figure). 

Finally we argue that the ideas may also be applicable to gauge theories in the form of world-sheet 
methods. Clearly the subject is still in its infancy and it would be exciting if some of the ideas can 
lead to important computational break-through. 
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